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Density-functional calculations are presented for high- 
pressure structural phases of S and Se. The structural phase 
diagrams, phonon spectra, electron-phonon coupling, and su- 
perconducting properties of the isovalent elements are com- 
pared. We find that with increasing pressure, Se adopts a 
sequence of ever more closely packed structures (/3-Po, bcc, 
fee), while S favors more open structures (/3-Po, simple cu- 
bic, bcc). These differences are shown to be attributable to 
differences in the S and Se core states. All the compressed 
phases of S and Se considered are calculated to have weak to 
moderate electron-phonon coupling strengths consistent with 
superconducting transition temperatures in the range of 1 to 
20 K. Our results compare well with experimental data on 
the /3-Po— >bcc transition pressure in Se and on the super- 
conducting transition temperature in /3-Po S. Further experi- 
ments are suggested to search for the other structural phases 
predicted at higher pressures and to test theoretical results on 
the electron-phonon interaction and superconducting proper- 
ties. 
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I. INTRODUCTION 

The chalcogens exist in a rich variety of crystalpStruc- 
tures at ambient pressures and upon compression, uu At 
low pressures, S, Se, and Te crystallize in different se- 
quences of insulating two-fold-coordinated molecular or 
polymeric structures. With enough compression, how- 
ever, all three elements transform into the same body- 
centered orthorhombic (bco) crystal structure. The bco 
structure is comprised of puckered layers with each site 
having a coordination number of four. With further com- 
pression, the more three-dimensional six-fold coordinated 
rhombohedral /3-Po structure becomes stable in all three 
systems. At even higher pressures, Se and Te are ob- 
served to transform from the /3-Po structure to the bcc 
structure, continuing the trend of increasing coordination 
and closer packing with pressure. Though is has been 
speculated that S should also undergo a /3-Po — > bcc tran- 
sition, recent first-principles calculations predict /3-Po S 
to become unstable .-first to the relatively open simple- 
cubic (sc) structured The calculations find simple-cubic 
S to be favored over a large range of pressures, with the 



bcc phase eventually becoming stable at very high pres- 
sures. 

Electrical measurements show that S, Sc, and Te are 
all metallic in the bco structure. Further, the metal- 
lic bco phases of S, Se, and Te are observed to be su- 
perconducting with transition temperatures of approx- 
imately 10, 5, and 3 K, respectively.Eltl Recent mag- 
netic measurements show that the superconducting tran- 
sition temperature in S jumps abruptly to kZ-iK upon 
transformation to the /3-Po structural phaseJij This is 
among the highest transition temperatures observed in 
elemental solids. First-principles calculations indicate 
that this high transition temperature results from a mod- 
erate electron-phonon coupling combined with a large av- 
erage phonon frequency, which sets the scale for the tran- 
sition temper ature.tl Earlier calculations had predicted 
the hypothetical bcc phase of S to be superconducting 
with a similar transition temperature of about 15 K near 
550 GPa.y As yet, no measurements of Tc have been 
reported on the higher-pressure phases of the heavier 
chalcogens. 

In this paper, we report density-functional calculations 
of the structural and superconducting transitions in S 
and Se in the /3-Po and higher- pressure phases. Though 
the two materials are chemically similar, they display 
striking differences in their properties under pressure. 
In addition to the simple-cubic structure's viability as 
a high-pressure phase for S but not Se, the stability of 
the fee structure is another difference between the ma- 
terials. An analysis of the contributions to the total en- 
ergy shows that these structural differences can be un- 
derstood in terms of differences between the S and Se 
cores. Although some of the trends in the calculated 
electron-phonon coupling and superconducting proper- 
ties are similar in the two materials, the superconduct- 
ing transition temperature is higher in S than in Se, even 
within the same structural phase. A comparison of the 
phonon spectra and electron-phonon coupling in these 
materials helps to explain the differences. 

The remainder of this paper is organized as follows. In 
Sec. II, we compare and contrast the calculated high- 
pressure structural phase diagrams for S and Se. In 
Sec. Ill, the calculated vibrational spectra and electron- 
phonon coupling parameters for the high-pressure phases 
are presented and discussed. We also discuss ways in 
which our results can be tested and compared with 
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diamond- anvil-cell-based experiments, including mea- 
surements of superconducting properties, transport co- 
efficients, and optical conductivity. Concluding remarks 
are given in Sec. IV. 

II. HIGH-PRESSURE STRUCTURAL PHASES 
A. Method 

Zero-temperature structural energetics are calculated 
using the plane-wave pseudopotential method within the 
local density approximation. The pseudopotentials are 
generated using the Troullier-Martins method.E3 Kinetic- 
energy cutoffs of 70 Ry and 40 Ry are used for the 
plane-wave expansion of the Kohn-Sham orbitals in S 
and Se, respectively. The correlation potential is ap- 
proximated using the Peijdew-Zunger parameterization 
of the electron-gas resultsjLj and the nonlinearity of the 
exchange and correlation interaction between the core 
and .valence charge densities is treated using a partial 
core.E-3 The Brillouin zone for each structure is sampled 
on Monkhorst-Pack meshes of at least 20 3 points Jl3 

B. Results and discussion 
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FIG. 1. Total energy versus rhombohedral angle a r at 
constant volumes for S and Se. At moderate compression 
both elements favor the /3-Po structure (a r w 104°); increased 
compression results in the sequence /3-Po^sc^bcc for S and 
/3-Po^bcc— >fcc for Se. 



The high-pressure crystal structures considered here 
for the chalcogens can all be described in terms of a 
single-atom rhombohedral unit cell with rhombohedral 
angles a r of 60° (fee), -00° (sc) , approximately 104° (/3- 
Po), and 109.47° (bcc)£3 Figure [l] shows the total energy 
for S and Se calculated as a function of the rhombohedral 
angle at constant unit-cell volumes. The curves for the 
two materials arc strikingly different, leading to differ- 
ent sequences of high-pressure phases. In particular, for 
the range of volumes shown in Fig. [I], the simple-cubic 
structure is never competitive in energy with the other 
structures considered for Se, while the fee structure is 
never favorable for S. 

For the largest unit-cell volume shown, correspond- 
ing to roughly 80 GPa in Se, the /3-Po structure with 
a r w 104° is lowest in energy. At a smaller volume of 
76 a.u. (roughly 140 GPa), the local minimum for /3-Po 
Se has disappeared, and the bec structure is favored. At 
even smaller volumes, the minimum at a r = 60° drops 
below that at 109.5°, indicating a bec to fee transition. 
Comparison of the enthalpies H = E to ta\ +pV calculated 
for the different structures yields a transition pressure of 
120 GPa for the /3-Po— >bcc transition, which is in rea- 
sonable agreement with the experimental value of 140 
GPa. For the bec— >fcc transition, we predict a transition 
pressure of about 260 GPa. 

For S, the /3-Po structure, which is favored at a unit-cell 
volume of 56 a.u., is calculated to become unstable under 
compression with respect to the simple-cubic structure, 



and the simple-cubic structure subsequently becomes un- 
stable to the bec structure (not shown in Fig. [I]). Us- 
ing enthalpies calculated for the static lattices, we find 
transition pressures of 260 and 540 GPa for these struc- 
tural transformations. Since we find that the structural 
phases of S differ significantly in their average phonon 
frequencies (see Sec. Ill), it is important to take into ac- 
count zero-point contributions to the energies. Inclusion 
of zero-point energies shifts the estimated pressures in S 
to about 280 GPa for the /3-Po— >sc transition and about 
500 GPa for the sc— >bcc transition. 

In terms of coordination number and packing, Se fol- 
lows the expected trend of increasing coordination with 
pressure, going from six (/3-Po) to eight (bec) to twelve 
(fee) . With each transformation, Se adopts a more closely 
packed structure. Sulfur, on the other hand, transforms 
from the six-fold coordinated (3-Po structure to the sc 
structure, which is also six-fold coordinated, but is more 
open. (The ratio of next-nearest to nearest-neighbor dis- 
tances is 1.24 and 1.41 in the /3-Po and sc structures, 
respectively.) While the more highly coordinated bec 
structure is adopted at large compressions, the even more 
densely packed fee structure is never favored. The key to 
these differences in the high-pressure phase diagrams of 
these isovalent elements lies in their cores. 

Figure || compares the relative stability of the sc and 
bec structures of S and Se as a function of volume. 
To facilitate the comparison, the volumes are measured 
with respect to the calculated equilibrium bec volumes 
of V$ cc = 98.6 and 125.9 a.u./atom for S and Se, re- 
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FIG. 2. Total energy difference between the simple cubic 
and bcc structures as a function of volume for S and Se. Also 
shown are the energy differences in the one-electron, Hartree, 
exchange-correlation (xc), and ion- ion contributions to the to- 
tal energy. Volumes are measured with respect to the theoret- 



ical zero-pressure volumes of V ' 



i.6 and 125.9 a.u./atom 



calculated for bcc S and Se respectively. 

spectively. The total energy can be divided into the 
one-electron contribution, which arises from the nonin- 
teracting kinetic energy and the electron-ion interaction, 
the electron-electron Coulomb energy, the exchange and 
correlation energies, and the ion-ion energy, which in- 
cludes the Ewald energy and a term that accounts for 
the difference between the pseudopotential and the pure 
Coulomb potential of the ions. In both materials, the re- 
pulsive Hartree and ion-ion terms favor the more closely 
packed bcc structure, which has a more uniform distribu- 
tion of both ionic and electronic charge. The attractive 
electron-ion Coulomb interaction favors the more open 
sc structure which tends to have very nonuniform charge 
distributions. The difference in total energy between the 
two structures is smaller than differences in the individ- 
ual contributions to the total energy. Figure || shows that 
the one-electron contribution in S favors the sc structure 
significantly more strongly than it does in Se, thereby 
tipping the balance to stabilize the open sc structure. 
Further decomposition of the one-electron contribution 
shows that the kinetic energy favors the more uniform 
bcc structure. It is thus the electron-ion interaction that 
stabilizes the sc S phase. 

The importance of the electron-ion interaction in pro- 
ducing different stable structures in S and Se can be un- 
derstood in terms of the difference in their cores. In 
particular, since the Ne core of S contains only s and p 
electrons, the 3d states in S have no orthogonality con- 
straint with the core, resulting in a strongly attractive d 
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FIG. 3. Electronic density of states of bcc and sc S cal- 
culated at a volume of 50 a.u./atom. The Fermi level is at 
eV. The d-projected density of states, plotted in lighter lines, 
shows that the sc structure has broader d-bands and a higher 
occupation of d states. 



pseudopotential in the core region. On the other hand, 
in Se, which has d states in the core, the repulsive Pauli 
component of the pseudopotential largely cancels the at- 
tractive Coulomb component in the core region, resulting 
in a relatively weak d potential. 

Under ambient conditions, the valence d states are un- 
occupied in both S and Se, but with compression, the d 
bands broaden and eventually cross the Fermi level. The 
electronic density of states for bcc and sc S are shown in 
Fig. [| The volume of 50 a.u./atom corresponds to about 
250 GPa, which is close to the pressure at which the sc 
phase is calculated to become stable. The thin lines show 
the ci-projected density of states. Since the d bandwidth 
varies roughly as C/d„ n , where C is the coordination 
number and d nn is the nearest neighbor distance,E3 the 
shorter bonds in the sc structure more than compensate 
for the lower coordination, leading to broader d bands in 
the sc phase as compared to the bcc phase. With more 
s-d transfer of electrons in the sc phase than in the bcc 
(or fee) structures, the deep S d pseudopotential becomes 
important in the energetics, stabilizing the sc phase. In 
Se, the sc phase similarly has larger d occupation than 
the fee and bcc phases, but because the d potential is 
weak, this s-d transfer does not lower the one-electron 
energy enough to stabilize the sc structure. 

Similar reasoning helps explain the absence of the 
close-packed fee structure in high-pressure S. Fig. ^ com- 
pares the relative stability of the fee and bcc phases of 
S and Se as a function of volume. Again the different 
structural preferences of S and Se can be attributed to 
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FIG. 4. Total energy difference between the fee and bee 
structures as a function of volume for S and Se. Also shown 
are the energy differences in the one-electron, Hartree, ex- 
change-correlation (xc), and ion-ion contributions to the total 
energy. Volumes are measured with respect to the theoreti- 
cal zero-pressure volumes of Vo = 98.6 and 125.9 a.u./atom 
calculated for bec S and Se respectively. 



III. PHONON SPECTRA, ELECTRON-PHONON 
COUPLING, AND SUPERCONDUCTIVITY 

A. Method 

The density-functional linear response methodEl is 
used to compute the phonon spectra for each struc- 
ture at several different pressures. The dynamical ma- 
trices are computed using the scheme described in Ref. 
pp| , where the self-consistent change in the Hamiltonian 
caused by ionic displacements is obtained by solving a 
Bethe-Saltpeter equation for the change in the charge 
density. The dynamical matrix is computed on a grid of 
6 3 phonon wavevectors q for (3-Po S and 8 3 wavevectors 
for sc and bec S. For Se, 8 3 , 10 3 , and 12 3 q points are 
sampled in the full /3-Po, bec, and fee Brillouin zones, 
respectively. 

The matrix element for scattering of an electron from 
state nk to state n'k' by a phonon with frequency w qj/ 
and eigenvector e qi/ is given by 



<?(nk, 7i'k', qy) = 



2Mw„ 



ik|eq„ • VRVscFlnV), (1) 



where VrVscf is the gradient of the self-consistent po- 
tential with respect to atomic displacements. This scat- 
tering gives rise to a finite phonon lincwidth, 



differences in the one-electron contributions to the total 
energy. With increasing compression, the one-electron 
energy in S increasingly favors the bec structure, while 
in Se it increasingly favors the fee structure. In both 
materials, the kinetic energy favors the more uniformly 
distributed fee structure while the electron-ion term fa- 
vors the bec structure. With the deeper d pseudopoten- 
tial in S, the importance of the electron- ion term in S 
is enhanced, stabilizing the bec structure. The differ- 
ence in openness between the fee and bec structures is of 
course much smaller than that between the sc and bec 
structures, so in comparisons of the fee and bec struc- 
tures, the energy differences are smaller and the energy 
balance is more subtle. 

This discussion of the differences between S and Se 
and the important role of the presence of d states in the 
core suggests that the heavier chalcogen Te should be 
similar to Se in terms of high-pressure structural phases. 
Indeed, like Se, Te is observed to transform from the j3- 
Po structure to the bec structure with pressured In our. 
calculations, as well as earlier ones by Kirchhoff et al.,tij 
this transition in Te is found close to the experimentally 
measured transition pressure of 27 GPa. In addition, we 
find that the close-packed fee structure becomes stable 
above about 80 GPa. Finally, as in Se, the sc structure 
is calculated to be energetically unfavorable as a high- 
pressure phase of Te.EI 



7qt , = 27r Wqi 4V(£ F )] 2 (<|.g ql ,| 2 )), 



(2) 



where N(Ep) is the electronic density of states per spin 
at the Fermi level. The double brackets ((...}} denote a 
doubly constrained Fermi surface average as defined in 
Ref. £y. This scattering process also contributes to the 
effective mass of the electrons via the mass enhancement 
parameter A, obtained from the wavevector- and branch- 
dependent contributions 
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hn N(E f )uj, 



(3) 



q'' 



by summing over branches v and averaging over wavevec- 
tors q. 

The electron-phonon spectral function, which mea- 
sures the effectiveness of phonons of a given energy to 
scatter electrons on the Fermi surface, plays a central 
role in the Eliashberg strong-coupling theory of super- 
conductivity. The spectral function is given by 



2 FH = 



1 



2nN(E F ) ^ 



S(CJ- W q!/ ) 



7q^ 



(4) 



Within this framework, the mass enhancement parameter 
A is proportional to the inverse-frequency moment of the 
spectral function. 

In the density-functional linear-response method, 
VrVscf is computed in the process of determining the 
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TABLE I. Calculated values of structural, vibrational, and 
electron-phonon parameters for /3-Po, sc, and bcc S at various 
pressures. The electronic density of states at the Fermi level 
N(Ep) is in units of states/eV/spin/cell. 





P 


V 


Or 


N(E F ) 


(w) 


(win) 


A 




(GPa) 


(a.u.) 






(meV) 


(meV) 






160 


57.87 


104.0° 


0.148 


46.4 


37.7 


0.76 


/3-Po 


200 


54.22 


104.1° 


0.143 


50.6 


40.2 


0.76 




280 


48.77 


104.4° 


0.140 


57.2 


42.3 


0.61 




280 


47.80 


90.0° 


0.159 


71.8 


48.1 


0.52 


sc 


320 


46.14 


90.0° 


0.146 


76.5 


52.8 


0.44 




550 


39.77 


90.0° 


0.118 


91.8 


57.8 


0.35 


bcc 


550 


38.81 


109.5° 


0.136 


67.0 


50.0 


0.70 



dynamical matrices. Therefore we calculate the electron- 
phonon matrix elements on the same grid of wavevectors 
used in the phonon calculations. For calculations of the 
phonon density of states and electron-phonon spectral 
function, which involve summing over q points through- 
out the Brillouin zone, a Fourier interpolation procedure 
is used to obtain the dynamical matrix and its dissipa- 
tive part on a denser mesh of q points.EHl The doubly 
constrained Fermi surface average of g is computed using 
Bloch functions on dense meshes of at least 30 3 k and k' 
points in the full Brillouin zone, with delta functions at 
the Fermi level replaced by Gaussians of width of order 
0.01 Ry chosen to reproduce the valuejjof N(Ep) obtained 
using the linear tetrahedron method.Eil 

B. Results and discussion 

Tables Q and |TJ list parameters that characterize 
the phonon spectra and electron-phonon couplings in 
the high-pressure phases of S and Se. Included are 
the average phonon frequency (u) and the logarith- 
mic A- weighted average phonon frequency (u)\ n ) = 
^PEqi/VM^l/Eq.V]- Some of tne values for 
N(Ef) and A listed in Table Q differ by a few percent 
from those published in our earlier papeiu on S because 
we use the more accurate tetrahedron method for calcu- 
lating the electronic density of states in the present work. 

Within each structural phase, compression stiffens the 
lattice and weakens the electron-phonon interaction, as 
indicated by the rise in (ui) and fall in A. The pressure- 
dependence of A varies from structure to structure. In sc 
S, A varies approximately as l/(w) 2 , as might be expected 
from Eq. (3). In /3-Po S, /3-Po Se, and bcc Se, however, 
the pressure-dependence of A is much weaker. These dif- 
ferences can be understood by examining contributions 
to A from different phonon modes. 

In some of the structural phases, it is possible to iden- 
tify specific phonon modes that contribute strongly to A. 
In the /3-Po phase of Se, for example, there are strong 
anomalies in the phonon dispersion curves along the T 
to F direction (i.e., the (110) direction in terms of the 



TABLE II. Calculated values of structural, vibrational, 
and electron-phonon parameters for /3-Po, bcc, and fee Se 
at various pressures. The electronic density of states at the 
Fermi level N(Ef) is in units of states/eV/spin/cell. 
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ot r 


N(E F ) 
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(Win) 


A 




(GPa) 


(a.u.) 






(meV) 


(meV) 






60 


93.33 


104.1° 


0.180 


23.5 


18.8 


0.58 


/3-Po 


80 


87.37 


104.8° 


0.176 


25.5 


19.5 


0.54 




100 


83.20 


105.2° 


0.174 


27.4 


21.6 


0.50 




120 


78.73 


109.5° 


0.201 


26.6 


17.6 


0.89 




140 


75.80 


109.5° 


0.193 


28.3 


18.8 


0.76 


bcc 


160 


73.32 


109.5° 


0.187 


29.4 


17.6 


0.77 




220 


67.36 


109.5° 


0.182 


33.3 


21.2 


0.73 




260 


64.25 


109.5° 


0.175 


35.6 


21.4 


0.71 


fee 


260 


63.82 


60.0° 


0.183 


34.93 


25.23 


0.59 
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FIG. 5. Phonon dispersion curves in /3-Po Se at 60 GPa. 
Directions are specified in terms of reciprocal-lattice basis vec- 
tors. Soft phonon anomalies related to nesting of the Fermi 
surface are evident along the 1 to F line. 

reciprocal lattice vectors bi), as shown in Fig. ||. Similar 
but less pronounced anomalies exist in the /3-Po phase 
of S. These anomalies are associated with nesting of the 
Fermi surface by the wavevector q w 0.7(bi — b2), as 
confirmed by calculations of the geometric nesting factor 
((<5(k — k' — q))). The strong Fermi surface nesting af- 
fects not only the phonon frequencies, but the linewidths 
as well. The combination of low frequencies and en- 
hanced linewidths results in a large contribution to A 
from this part of the /3-Po Brillouin zone. Anomalies are 
also present in the phonon spectra of bcc S and bcc Se. 
Figure |^ shows the dispersion curves for bcc Se at 120 
GPa. While not correlated with nesting vectors, the low- 
frequency anomalies in the transverse branches along the 
F to N and T to H lines nevertheless give large contribu- 
tions to A. 

The soft phonon anomalies in the bcc and /3-Po phases 
are more weakly pressure-dependent than the overall 
phonon spectrum. The pressure dependence of A in these 
phases is strongly influenced by these persistent anoma- 
lies and is therefore weaker than the {lo)~ 2 dependence 
suggested by Eq. (3). In sc S, the electron-phonon cou- 
pling is not found to be dominated by any particular 
phonon mode, so the pressure dependence of A is to a 
good approximation determined by that of the average 
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FIG. 6. Phonon dispersion curves in bcc Se at 120 GPa. 
Directions and polarizations are specified in terms of recipro- 
cal-lattice basis vectors. Phonon anomalies are present in the 
(lll)-polarized branch along P to N, and in the transverse 
branches along T to H. 
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FIG. 7. Electron-phonon spectral functions calculated for 
/3-Po, sc, and bcc S. 



phonon frequency. 

As can be seen in Tables | and B transformations 
from structure to structure are usually accompanied by 
abrupt changes in the phonon and electron-phonon pa- 
rameters. In comparing the various structural phases of 
S, the simple-cubic phase stands out in having the weak- 
est electron-phonon coupling parameters. The spectral 
function for sc S, plotted in Fig.^, is strikingly different 
from those of bcc and /3-Po S. Both the average and max- 
imum phonon frequencies in the sc phase are much larger 
than those in the other phases. The higher frequencies in 
the sc phase are due to the shorter nearest-neighbor dis- 
tances and the strong concentration of charge along the 
bonds in the open sc lattice. With these stiff bonds, the 
(llO)-polarized transverse branch along the (110) direc- 
tion and both of the transverse branches along the (111) 
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FIG. 8. Electron-phonon spectral functions calculated for 
/3-Po, bcc, and fee Se. 

direction are nearly as high in energy as the longitudi- 
nal branches. This leads to a phonon density of states 
and a spectral function that are dominated by a high- 
frequency peak. The large phonon energy scale and the 
lack of low-energy modes that couple strongly to elec- 
trons, as evidenced by the relatively large value of (co\ n ), 
combine to drive A down in the sc phase of S. 

Among the high-pressure phases of Se, the bcc phase 
stands out in that A is largest in this phase. Although 
the average phonon frequencies are similar to those in the 
neighboring phases near the transition pressures, (uj\ n } is 
significantly depressed in the bcc phase. This is a sign of 
the importance of the low-frequency anomalies already 
identified in the bcc Se phonon spectra. These anomalies 
enhance the spectral function at low energies, as shown 
in Fig.||, and they have a significant effect on A because 
they appear at low frequencies. Although /3-Po Se also 
has soft phonon anomalies, there are fewer of them due to 
the lower symmetry of the lattice, and they have a smaller 
effect on A. No comparable soft phonon anomalies are 
found in the fee phase. 

In comparing analogous structural phases of S and Se, 
the differences in electron-phonon coupling strength re- 
flect not only differences in the phonon spectra, but also 
differences in the electron-phonon matrix elements. The 
stronger potential in S gives rise to a larger change in 
the bare potential VR,I4>are when atoms are displaced. 
Furthermore, because the d electrons are more tightly 
bound in S, these electrons are less effective in screening 
VR,Vbarc- In combination, these two effects result in sig- 
nificantly larger matrix elements of VrVscf in S than in 
Se. Combined with the smaller mass in S, this leads to 
larger linewidths. This is evident in both the /3-Po and 
bcc phases, where the phonon frequencies in S are a fac- 
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tor of two or more higher than those in Se, yet the mass 
enhancement parameters in the two materials are close. 

Our results can be compared to earlier density- 
functional-based studies of the electrcp.pl|ionon interac- 
tion in /3-Po and bcc Te and bcc S.EIrEJ Soft phonon 
anomalies very similar to those discussed here have been 
identified to be important pin the electron-phonon cou- 
pling in /3-Po and bcc TeO In addition, the electron- 
phonon mass enhancement parameter has been calcu- 
lated to undergo a large jump at the /3-Po to bcc tran- 
sition in Te, just as we find for Se. For bcc S, our value 
of A = 0.70 at 550 GPa is somewhat larger than the 
value of 0.58 reported in Ref. |ll| for a slightly higher 
pressure of 584 GPa. Since values of A qz , are in rea- 
sonable agreement along high-symmetry directions, it is 
likely that the difference arises from a difference in sam- 
pling of wavevectors: In the present work, a uniform grid 
of points throughout the Brillouin zone is sampled, while 
in Ref. [H], A q „ is calculated along a few high-symmetry 
directions and spherically averaged to estimate A. 



C. Making contact with experiments 

The need to carry out experiments inside high-pressure 
cells limits the types of experimental probes available to 
investigate the electron-phonon interaction in the metal- 
lic phases of the chalcogens. In particular, the most 
common probe of a 2 F in superconductors, quasiparticle 
tunncling,E£l is challenging because of the need to make 
well characterized tunnel junctiptis connected to leads 
that enter the high-pressure cellES 

Measurement of superconducting properties such as 
the transition temperature Tc, the thermodynamic crit- 
ical field He, the zero-temperature gap Ao, and the 
isotope effect on Tc, would be pae avenue for probing 
the electron-phonon interaction.EH Since superconduct- 
ing properties depend not only on the electron-phonon 
spectral function a 2 F, but also on the Coulomb pseu- 
dopotential fi* , it is necessary to measure a combination 
of superconducting properties to get a handle on both 
the electron-electron and the electron-phonon parame- 
ters. In the weak coupling regime, Tc, Ao, and He all 
essentially depend on the difference between A and p* , 
rather than on a 2 F and p* separately. In this case, Tc 
and the isotope effect on Tc, with their different a 2 F and 
p* dependences, would be the most useful combination 
of measurements. 

Figure 9 shows our calculated superconducting transi- 
tion temperatures for the high-pressure phases of S and 
Se for different values of p* , obtained by using the calcu- 
lated a 2 F functions as input to the Eliashbcrg equations. 
AvailabLe,experimental data for bco Sejj bco S,E!J and /3- 
Po S,Ejc3 are included. The theoretical data points in 
Fig. |^ correspond to /i*(w m ax) = 0.11, where w max is the 
maximum phonon frequency. This representative value of 
M*(w m ax) is chosen because it reproduces the experimen- 
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FIG. 9. Superconducting transition temperature of S and 
Se as a function of pressure. Squares, circles, up triangles, 
and down triangles represent results of calculations for the 
/3-Po, bcc, fee, and sc structures, respectively. The error bars 
show the range of Tc for 0.09 < /i*(w max ) < 0.13 and the 
data points correspond to /i*(oj max ) = 0.11. Experimental 
data are plotted as crosses. The dashed lines connecting cal- 
culated points serve as guides to the eye, and the vertical gray 
bars roughly separate the pressure ranges over which different 
structures are stable. 



tally measured T c of 17 K at 160 GPa in S.EJ The error 
bars show the effect of varying /x*(w max ) by ±0.02. While 
there is no reason to assume p* remains constant with 
changes in pressure or structure, or that p* in Se should 
be the same as in S, p* typically lies within the range of 
0.1 to 0.14 for a wide range of materials. Furthermore, 
a rough upper bound for p* given by [ln(_Ep /w ma x)] _1 
suggests p* is less than about 0.16 in all the phases con- 
sidered here. 

The calculated transition temperatures follow the same 
trends as A, with large increases in Tc upon transition 
to both bcc Se and bcc S, and downward jumps in Tc 
at the bcc^fcc transition in Se-,and the /3-Po^sc tran- 
sition in S. As reported earlier,Q our analysis attributes 
the large Tc of 17 K in /3-Po S not so much to strong 
electron-phonon coupling, but rather to a combination 
of moderate coupling and a large phonon energy scale. 
For analogous structural phases in S and Se, the primary 
reason for the higher transition temperatures in S is the 
stiffcr lattice, which sets the energy scale for Tc- 

The high phonon frequencies suggest that vertex cor- 
rections might be important in these materials. The 
other determining factor for the strength of vertex 
corrections _is the magnitude of the electron-phonon 
interaction,E3 which is moderate in both S and Se. In- 
deed, we estimate that the effect of the vertex corrections 
on Tc is on the order of one percent, and can be neglected 
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FIG. 10. Normal-state (a) reflectance and (b) optical con- 
ductivity calculated for /3-Po S at 160 GPa and T=0 K. Re- 
sults are shown for different impurity scattering times r. In 
(c), the optical conductivity with r = 0.1 eV is redrawn along 
with the corresponding Drude conductivity (tq(uj) assuming 
scattering by impurities only. Also plotted is cro[w(l + A)], 
showing that electron-phonon scattering narrows the Drude 
peak by a factor of 1 + A. 



given the uncertainty in /i*. 

Although we have used the Eliashberg formalism to 
calculate Tc for all the phases, the coupling is sufficiently 
weak that the Allen-Dynes approximate formula for TcEj 
gives similar results. Furthermore, the ratios 2Ao/ksTc 
and jT^/H^, where 7 is the linear coefficient of the spe- 
cific heat, are calculated to be close to the BCS values of 
3.53 and 1.68, respectively. For /3-Po S at 160 GPa, for 
example, we find 2A /k B T c = 3.70 and jT c /H c =1.66, 
assuming /x*(w max ) = 0.11. 

As suggested above, measurement of the isotope effect 
on Tc would be one way to obtain more information on 
the Coulomb parameter. For S at 160 GPa, we estimate 
that the isotope exponent a = —d\nTc/d\nM changes 
from 0.48 to 0.43 as /x* is varied from 0.09 to 0.13, as- 
suming isotopic masses M of 32 and 36. To distinguish 
between /i* = 0.09 and /1* = 0.13 then requires being 
able to measure Tc to an accuracy of better than 0.1K. 

Infrared spectroscopy can be used to probe a 2 F in- 
dependent of ^*E3E3 In the normal state, the phonons 
contribute to the optical conductivity via their renormal- 
ization of the electronic quasiparticles and by phonon- 
assisted scattering processes. Panels (a) and (b) in 
Fig. |h] show the normal-state T — reflectance and 
optical conductivity calculated for different values of im- 
purity scattering rates for /3-Po S at 160 GPa. The ex- 



pected Holstein structure, which occurs at frequencies 
near the phonon frequencies, is small because the cou- 
pling strength is relatively weak. Hence, it is likely that 
the extraction of a 2 F from phonon structure in the op- 
tical conductivity data would be challenging for this sys- 
tem. However, optical conductivity data could still be 
useful for estimating A. The electron-phonon interaction 
narrows the Drude peak associated with impurity scat- 
tering by a factor of 1+A due to the renormalization of 
electronic quasiparticle energies and weights and trans- 
fers the low-frequency weight to the mid-infrared region, 
as illustrated in Fig. |l0|(c). With an independent de- 
termination of the plasma frequency u> p , the width and 
height of the Drude peak could be used to estimate A. 

Perhaps a simpler way to determine A from experi- 
ments is via the dc conductivity. At temperatures on 
the order of the Debye temperature, the phonon-limited 
electrical resistivity is linear in temperature with a slope 
proportional to A tr = 2 J duja 2 r F(uj)/uj. Like a 2 F, the 
transport spectral function a 2 r F measures the effective- 
ness of phonons to scatter electrons on the Fermi sur- 
face, but it is weighted to take intp-,account the change 
in direction of the electron velocity.til Experimental mea- 
surements of the electrical resistivity have been reporto 
for the bco and lower pressure phases of S and Se.B 
Extensions of these measurements to higher pressures 
would test our results for the electron-phonon coupling 
strength. Of particular interest is the linear coefficient 
of the resistivity in /3-Po S. We calculate Atr = 0.78 at 
160 GPa, close to our result for A. We believe this is the 
most direct way to verify our theoretical prediction that, 
despite the large Tc, the electron-phonon coupling is not 
particularly strong in this phase. 



IV. CONCLUSIONS 

We have presented density-functional based calcula- 
tions on the pressure-induced structural phase transitions 
in S and Se. Although there are some similarities in 
the high-pressure phase diagrams of these isovalent ele- 
ments, there are also striking differences. With increasing 
compression, Se adopts a sequence of ever more closely 
packed structures (/3-Po^bcc^fcc), while S favors more 
open structures (/3-Po^sc^bcc). These differences can 
be understood in terms of the deeper d pseudopotential 
in S arising from the lack of d states in the S core. 

All the high-pressure phases of S and Se are calculated 
to have moderate electron-phonon coupling strengths and 
superconducting transition temperatures in the range of 
0.5-20 K. Structure and pressure dependences of the cou- 
pling strength and Tc can be understood in terms of 
changes in the phonon spectra. In particular, the large 
observed Tc in /3-Po S is calculated to arise from a com- 
bination of a moderate mass enhancement parameter A 
and a large overall phonon energy scale. The large in- 
crease in A and Tc predicted upon transformation to the 
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bcc phase in both S and Se is due to the presence of soft 
modes in the bcc phonon spectra that strongly couple to 
electrons. 

We suggest that measurements of the electrical resis- 
tivity as a function of temperature, the transition tem- 
perature Tc, and the isotope effect on Tc would be the 
most promising avenues for testing our results on the 
electron-phonon interaction in these compressed phases. 
Further, Tc measurements could be used to search for the 
predicted sc phase of S and the predicted fee phase of Se 
since our results indicate that Tq should change abruptly 
when these structural transformations take place. The 
calculated transition pressures for both sc S and fee Se 
are within the range of current diamond-anvil-cell exper- 
iments. 
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